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The possibility that the Galactic dark matter is composed of neutralinos that are 
just above half the Z° mass is examined, in the context of the Galactic positron 
excess. In particular, we check if the anomalous bump in the cosmic ray positron 
to electron ratio at 10 GeV can be explained with the "decay" of virtual Z° 
, bosons produced when the neutralinos annihilate. We find that the low energy 

O behaviour of our prediction fits well the existing data. Assuming the neutralinos 

I annihilate primarily in the distant density concentration in the Galaxy and allowing 

O ' combination of older, diffused positrons with young free-streaming ones, produces 

a fit which is not satisfactory on its own but is significantly better than the one 



obtained with homogeneous injection. 



1. Introduction 



■ The possibility that weakly interacting dark matter particles (WIMP's) 

could annihilate into detectable cosmic radiation was suggested by Silk and 
Srcdnicki . Tylka and Eichler 2 noted a reported positron excess, curiously 
localized near 10 GeV 3 > 4 ' 5 , and considered whether it could be due to the 
annihilation of photinos (as a simple example of neutralinos) in the tens of 
GeV mass range. The difficulty was that this process, given the laboratory 
constraints on the neutralinos, seemed to fall short of providing enough 
positrons, and the results were not published. Nevertheless, various papers 
on this excess eventually appeared, and noted that the potential for positron 
excess could be bolstered by dumpiness in the annihilating dark matter or 
by decay of weakly unstable dark matter particles. 

The approach usually found in present literature is to try and fit the 
overall e + /(e + + e~) ratio, without giving special attention to the curi- 
ous behaviour at 10 GeV 6 . Baltz and Edsjo 7 considered a whole class 
of minimal standard supersymmetric models and failed to get any non- 
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monotonicity in the e + /(e + + e~) ratio. They note that an enhancement 
factor, presumably due to clumping, of at least 30 to 100 is needed to 
account for the observed positron excess. 

Eichler and Maor 8 considered the possibility that the observed excess 
in positrons, and in particular the anomalous behaviour around ~ 10 GeV, 
is a result of annihilating dark matter particles with mass just above mz/2. 
A non-relativistic virtual Z° decay (i.e. when the rest mass of the annihi- 
lating dark matter particle is slightly above 1/2 the Z° mass) provides a 
remarkably good fit to the observed e + /(e + + e~) ratio below 10 GeV. At 
higher energies, however, the predicted e+/(e+ + e~) ratio rises above the 
observed values, within conventional assumptions about the injection and 
propagation. In particular, it was assumed that the positrons and primary 
electrons are each injected homogeneously, with the same spatial profile, 
and that their propagation in the Galaxy is identical. Relaxing the as- 
sumption of homogeneous injection for the Z° decay products was explored 
in a consequent work 9 , and the fit to the data at high energies improved 
considerably. A summary of the virtual Z° decay model is presented here. 

2. The basic equation and its solutions 

The fact that the observed positron excess appears at the characteristic 
energy range of positrons from non-relativistic Z° decay motivates our con- 
sideration of this origin. Such Z° particles could be produced by the anni- 
hilation of neutralinos that are just beyond half the Z° mass (so that they 
do not contribute to the Z° decay width). Alternatively, they could be 
produced via the slow decay of some other particle that is just higher than 
twice the Z° mass. One could even imagine that such a particle would cou- 
ple to non-relativistic matter but not to relativistic matter, e.g. a massive 
dilaton coupled to the trace of T^, so that it would decay preferentially 
into Z's, but not directly into the lighter decay products of the Z. 

In order to achieve enough annihilation of the neutralinos, dumpiness 
needs to be invoked. This is a generic problem of models which make use of 
dark matter annihilation products. The fact that clumping is required, and 
that a likely place for this is near the Galactic centre, brings us to consider 
that the positrons in our neighborhood that are dark matter annihilation 
products would have a minimum age. The time needed to diffuse from the 
source to our neighborhood, and the age distribution of the positrons that 
make it to the Earth's vicinity will result in fewer young positrons than 
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one associates with the standard leaky box model. We therefore consider, 
in addition to standard assumptions about propagation, that the positrons 
are injected by an effectively point source at a finite distance, and look at 
the diffused equation. 

The steady state diffused equation for the particle number density is 

^ = Q = bn-Rn+— A (mz^n) + I(x)6(r) (1) 
at m z ox \ at J 

x = E/mz, D is a diffusion operator, R = Bx 05 is the escape rate, mzj§~ 
Ax 2 with A = 8.5 x 10~ 16 erg/s is the Compton loss rate, corresponding 
to an electromagnetic energy density in the Galaxy of 10~ 12 erg /cm 3 and 
I{x) is the injected spectrum. 



For a homogeneous injection, Dn = 0, 
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A one-dimensional diffusion operator with coefficient T>, D = T>-^, and 
with boundary conditions such that ^| r=i = (conserving the number of 
particles except for the escape term), gives 
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While V (the diffusion coefficient) and L (the size of the leaky box) are 
free parameters, we took r = 8 Kpc, the distance to the galactic centre. 
K = D/r 2 gives the inverse time for diffusion. 



3. Injected spectrums 
3.1. Backgrounds 

The injected spectrum of primary electrons is taken to be I(x) = Cx~ 2 . 
For the positron background we take Dx~ 2 8 , presumed to come from cos- 
mic ray collisions. 
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3.2. Z° Products 

We now calculate the spectrum of positrons resulting from the annihilation 
of the neutralinos (x) via the virtual Z° channel assuming that m x is just 
above m^/2. It is essentially Z° decay. We have taken three generations of 
particle families, except for the top quark since it is much heavier then the 
mass of the Z° . All other particles were considered massless, and the cal- 
culation is in zeroth order. The positrons' as well as the electrons' injected 
spectrum has contributions from the following channels: 

• Z — > ee 

The branching ratio of direct decay to electron-positron pair is 



• Z — > fijl 

All the produced muons eventually decay into z/V e e, in an ex- 
change of a W boson. The resulting spectrum energy 10 in the 
muon rest frame is: 



where e is the positron's energy in the muon's rest frame. To boost 
the spectrum into the observer's frame, we have taken j3 ~ 1 and 
7 = m z l '(2m^). Assuming the muon decay is isotropic: 



The tau can decay into 3 colors of ud pairs, into electrons directly, 
or into muons which then decay into electrons. For the 20% of taus 
which decay directly into electrons, the calculation is the same as 
the above, with a branching ratio of 0.0344 x 0.2. 

For the 20% of the r — > a — > e channel, the resulting contribu- 




/(e) = £ (3m^ - e±) 





with a similar branching ratio, 0.0344. 



tion is: 




with a branching ratio of 0.0344 x 0.2. 
We have neglected the 60% of the t — > q. 
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• Z -> qq 

The hadronic channel of the Z decay eventually contributes some 
electrons and positrons through the production of pions. We have 
estimated the energy spectrum of the pion production from 11 : 

P{E V ) = I0 afc - bfc2 
oi = 3, bi = 10 < x < 0.1 

a 2 = 2, b 2 = 4 0.1 < x < 1 

with x = 2E n /mz- 

In the pion's rest frame, the energy spectrum of the electrons is 
close to a delta function, with an energy of ~ 45 MeV . In trans- 
ferring back to the lab frame, we have used (3 w 1 and 7 = E^/m^. 
The resulting contribution is: 




The branching ratio of of the Z° decay into hadrons is 0.6916. 

After taking the right weight of each of these channels, the injected 
spectrums are: 

I{x) = Qm\QNI h {x) + 0.03447V x (4) 

\( 1\ 70 168 , 6272 , /10 72 , 256 A, ,„ ." 

5 x H x 2 H x 6 1 x z H x A In (2.x) 

V 2 J 27 5 135 V 9 5 45 J K ' 

where N is the annihilation rate density. The needed injection rate to 
supply the observed positrons is about 1 x 10 _29 cm _3 s _1 above 10 GeV. 
The actual annihilation rate density N therefore needs to be about (2.2 x 
0.034)" 1 times that, N = 1.3 x 10" 29 cm -3 s _1 . Here the factor of 2.2 
assumes that Z decay into electrons, muons, and 20 percent of the taus 
results in positrons above 10 GeV. The need for significant dark matter 
clumping has not been avoided in the present scenario. 

4. Results 

We first consider a homogeneous injection of both backgrounds and the Z° 
products, see Figure 1. A striking feature of this figure is the good fit below 
10 GeV, which is mainly due to the positrons that emerge from decaying 
muons. Above 10 GeV, however, the fit grows very quickly, which is not 
compatible with the plateau that the data suggests. The reason for this 
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rise is that some Z 0, s decay directly into high energy e + e~ pairs so that 
the e + energy is half the Z° mass, and this gives rise to a high energy bump 
in the e + /(e + + e~) ratio at about 50 GeV. Even though the branching 
ratio into direct e+e~ is rather small (0.0344), it is dominant over the 
power-law background contribution at this energy. While this bump can 
be partially washed out by losses and escape, it was found that the high 
energy e + /(e+ + e~) ratio is nevertheless apparently too high to fit the 
observations. As discussed in Eichler (1989) this is a generic problem for 
any positron source that is significantly harder than the primary electrons 
above 10 GeV. 

Next we consider an inhomogeneous injection for the Z° products, while 
still keeping the backgrounds homogeneous. Figure 2 shows a fit with a 
single point source of Z° decay and a single diffusion coefficient. The good 
fit to the low energies from Figure 1 is still present, but the excess in energies 
toward x = 1/2 is now is suppressed by the finite age effect. As the figure 
shows, we are now facing a scenario which is opposite to the homogeneous 
injection case; the finite age effect tends to suppress the high energy excess 
at the price of killing it off altogether. 

However, there are several possibilities that dispel this problem: There 
may be more than one source, and there may be more than one route by 
which the particles diffuse or free-stream from the source to our vicinity. 
High energy particles diffuse much less than low energy ones because they 
are fewer in number and create less waves. So their self-generated scattering 
is less efficient. Thus, the fraction of free streaming particles should be 
higher at higher energy. Figure 3 shows a combination of two Z° decay 
components, an older, larger one that arrives via diffusion, and a younger, 
smaller component that has managed more free streaming. This figure 
illustrates that if one takes an age distribution into account, the flexibility 
in adjusting the high energy spectrum becomes much larger, and can be 
fitted to the data. 

For sake of comparison, we also include a power law injected spectrum. 
Figure 4 shows various power laws, and Figure 5 shows two components with 
different ages. We find that as long as the injected power law is hard enough, 
one can produce a low energy dip. While we find that the goodness of fit 
provided by Z° decay and by power law are qualitatively the same, we have 
deliberately not quantified this with the standard statistical measures. The 
GeV dip which we are trying to address formally has less statistical weight 
compared with the overall shape of the data, but it is the more striking 
observational result and the one that seems hardest to explain. Trying 
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to get the best parameters (for either power law or Z° decay as injected 
spectrums) would wash out the low energy behavior that we are focusing 
on. As we are interested in exploring the possibility of connecting the GeV 
behavior to properties of dark matter more than constraining parameters, 
we settle for the time being in showing characteristic examples. 

5. Discussion 

We have considered a scenario for positron excess via non-relativistic Z° 
decay in which the energies of the decay products are as low as possible and 
the neutralino annihilation rate that would give rise to the Z 0, s as high as 
possible. We considered homogeneous injection, as well as injection from a 
point source at a finite distance. 

Homogeneous injection fits the data remarkably well for energies below 
10 GeV, but produces too many high energy positrons (originating from 
the direct Z° — > e+e~ channel). A single source with a single diffusion 
coefficient does not provide a good fit to the high energy either, and the 
most promising scenario is to allow a small fraction of the positrons to free 
stream and arrive at the Earth's vicinity much younger than the rest. 

Although we can reproduce the 7 GeV dip, the peak at E <~ 15 GeV is 
still too big for the HEAT data. This seems to be a generic feature of our 
results, regardless of whether the injection source is virtual Z° decay or a 
power law. The problem would be worse if the virtual Z° had an energy 
well above mz- Knowing that the virtual Z° must be close to its mass 
shell if it is to provide a decent fit suggests that its loop corrections would 
be large and it might be discernable or falsifiable with particle collider 
data. Further possible observational consequences will be considered in a 
forthcoming publication 9 . 
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Figure 1. / = e + /(e + + e~) as a function of x = E/mz- Our model compared 
with HEAT 94+95 data The data and its error bars are marked with the crosses. The 
parameters of this example arc A = 8.5 X 10 -16 erg/s, B = 7.6 X 10~ 15 1/s, C = 
9.1 X 10~ 28 e/(cm 3 ■ s), and D = 1.3 X 10~ 31 p/(cm 3 ■ s). 
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Figure 5. The e+/(e+ +e ) as a function of x = E/mz, for a combination of 2 sources 
of power law (w = -0.3) injected spectrum. A = 8.5 X 10~ 16 2I£, B = 4.4 X 10~ 15 ±, 
C = 1.7 x 10- 29 —k-, D = 1.1 x 10" 31 -A— , = 7.2 x 10~ 17 ± and K 2 = 
1.2 x 10 -16 -. The ratio between the two components is 10 : 1. 



